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Dynamic effects induced by renormalization in anisotropic pattern forming systems 
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The dynamics of patterns in large two-dimensional domains remains a challenge in non-equilibrium 
phenomena. Often it is addressed through mild extensions of one-dimensional equations. We show 
that full 2D generalizations of the latter can lead to unexpected dynamical behavior. As an example 
we consider the anisotropic Kuramoto-Sivashinsky equation, that is a generic model of anisotropic 
pattern forming systems and has been derived in different instances of thin film dynamics. A 
rotation of a ripple pattern by 90° occurs in the system evolution when nonlinearities are strongly 
suppressed along one direction. This effect originates in non-linear parameter renormalization at 
different rates in the two system dimensions, showing a dynamical interplay between scale invariance 
and wavelength selection. Potential experimental realizations of this phenomenon are identified. 

PACS numbers: 05.45.-a,47.54.-r,68.35.Ct 
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The self-organized formation of patterns in non- 
equilibrium systems is a fascinating topic that has fo- 
cused a large attention in the last decades. Examples 
range from galaxy formation to sandy dunes, to nanos- 
tructures [1]. While regular patterns like stripes and 
hexagons that are characterized by a single length-scale 
£ are well understood when the lateral system size L is 
comparable to ^, their dynamics becomes much more 
complex in the large domain limit L ^ i. Indeed, in- 
tricate structures ensue, like spatiotemporal chaos, spi- 
ral waves or quasiperiodic patterns [2]. Another source 
of complexity derives from dimensionality. While a uni- 
fied description of one-dimensional patterns is available 
through the Ginzburg-Landau equation [i| , this is not the 
case for two-dimensional systems, strongly anisotropic 
problems providing particularly challenging cases. 

Many systems depending of two space variables are 
studied within the context of 3D localized structures [3| , 
like vortices in plasmas [4], or solitary waves in fluids 
[5|, [6| . The dynamic equations considered are frequently 
mild extensions of ID equations in which only specific 
terms are turned into 2D operators. This allows to probe 
the behavior of a given localized structure when the spa- 
tial dimension is increased. In other cases, the equation 
derives from first principles, as when studying anisotropic 
surface tension and kinetics in solidification systems [7]. 
A prototype model appearing in all these studies is the 
anisotropic Kuramoto-Sivashinsky (aKS) equation 



dth = v^jcdlh ^ Vyd^h ^ -^ {dxhf 



^ {dyhf 



-K^dlh - 2}C:,ydtd;h - Kyd^h + r? 



(1) 



For the sake of definiteness, we will keep in mind a phys- 
ical picture in which /i(x, ?/, t) is interpreted as the height 



of a surface above point (x, y) on a reference plane at 
time t. Indeed, particular instances of Eq. ([1]) have been 
derived in various contexts of thin film dynamics like sur- 
face nanopatterning by ion-beam erosion H, epitaxial 
growth [9-11] or solidification from a melt (l2(. In Eq. 
([1]), the morphological instability leading to pattern for- 
mation is implemented by the coefficients ^x^y, at least 
one of them being negative. Terms with coefficients Kj 
provide dissipation at the smallest scales, while the non- 
linearities proportional to Xx^y stabilize the system. We 
have incorporated a Gaussian, zero-average, uncorrelated 
noise r]{x^y,t) as a means to explore the aKS equation 
in a large domain. Indeed, for L ^ i^ the determinis- 
tic KS equation is well known to display spatio-temporal 
chaos and a steady state with strong height fluctuations 
[l3l , \1^ . Introduction of noise helps to elucidate these for 
times after onset of the morphological instability |15l-[l7j , 
while it is not essential for the occurrence of the morpho- 
logical transition we are studying in this paper. 

Regarding Eq. ([1]) as a model of 3D localized struc- 
tures, e.g. the description of solitary fluid waves moving 
down an inclined plane corresponds to Uy = Xy = 0, 
JCx = ^y = ^xy [y, [13 5 while in the solidiflcation system 
i^y = —fx^ Xy = Ky = JCxy = [?] . Tfius Eq. ^ is 
tailored to preserving the quasi one-dimensional features 
of speciflc solutions. In this work, we show that small 
deviations from flne tuned conditions such as the latter 
are able to induce dynamic effects in the system that 
unavoidably require a full two-dimensional description. 
Speciflcally, a rippled pattern appearing at short times 
along one of the system directions rotates by 90° during 
the evolution, leading at longer times to ripples oriented 
in the perpendicular direction. This dynamic transition 
occurs as a result of the different rates at which fluctua- 
tions renormalize due to the inhomogeneous strengths of 



nonlinear effects along the two space dimensions. 

We start by exploring numerically the behavior of Eq. 
([T]). It is convenient to first bring it to dimensionless 
form. As the properties to be studied are not conditioned 
by the anisotropy of the dissipative terms, we will restrict 
ourselves to the case in which JCx = JCy = JCxy = JC. 
Then, defining rescaled coordinates f = —{\iyx\^x/^)t^ 
x' = (|z/^|//C)i/2x, h' = -{\x/2ux)h leads to 

dth = -dlh - a^d^h + {d^hf + ax {dyhf- V^h + ^, (2) 

where primes are dropped, ^ is a noise term with a 
rescaled variance, and the ratios a^ = fy/vx and a a = 
^y/^x control the linear and the nonlinear anisotropics, 
respectively. We perform numerical simulations of Eq. 
([2|) for flat initial conditions by using both a finite differ- 
ence scheme with periodic boundary conditions pj|] , and 
alternatively a pseudospectral scheme [20,l2l|. 

Fig. [1] shows the time evolution of the surface rough- 
ness W"^{t) = ((1/^2) ^^{hr{t) - h{t)Y), where brackets 
denote average over noise realizations and bar denotes 
space average. Two different parameter conditions are 
considered for ajy < 1 [22]. In the first one (I), the non- 
linear couplings are comparable in the two directions x 
and y {ax = 0.5), while the second one (II) is represen- 
tative of conditions in which Xy is strongly suppressed 
(0 < ax ^ 0.1). Values of time at which W{t) changes 
behavior significantly are marked by showing the cor- 
responding surface morphologies. For both conditions, 
ajy < 1 induces a linear instability leading to forma- 
tion of a ripple structure with crests oriented parallel 
to the y axis. This takes place at time to in Fig. [TJ 
at which W{t) grows exponentially. The corresponding 
morphologies are statistically indistinguishable for both 
conditions, so that a single common snapshot (Ii,IIi) is 
shown in Fig. [TJ For later times, this regime is followed 
by non-linear stabilization inducing at time ti a slower, 
power-law growth rate for W{t). Morphologically, this 
type of growth is characterized by a progressive blur- 
ring of the early time pattern (morphology I2) and the 
dominance of height fluctuations associated with kinetic 
roughening [13|, [23|. For condition II, noise effects also 
dominate for time ti < t < t2 (morphology II2). At time 
t = t2 there is a second rapid increase in W{t) that mor- 
phologically corresponds to the formation of a new ripple 
structure, but now with crests parallel to the x axis. The 
wave-length of this new pattern is larger than that of 
the initial one, see morphology II3. Again, nonlinear ef- 
fects stabilize growth in amplitude for times t > ^3, so 
that W{t) displays kinetic roughening properties similar 
to those occurring for condition I at long times [23]. 

In the simulations, a a is always positive, so that the 
known cancellation modes of the aKS equation [11] are 
not solutions of Eq. ([2]). The dynamic behavior seen so 
far can be interpreted using results for the noisy ID KS 
equation. Thus, as already argued for by Yakhot [2^] 
in the deterministic limit, the noisy KS (nKS) equation 
is known to undergo a renormalization process through 
which the negative coefficient of the second order linear 




FIG. 1: Surface roughness vs time for conditions I (circles) 
and II (diamonds). Simulations of Eq. ([2|) are performed with 
a^ = 0.1, D = IQ-^ At = 5 X 10-^ L = 512, Ax = 1, and 
ax = 0.5 (I), 0.025 (II). W{t) is averaged over 100 noise real- 
izations. Top views of the surface morphology (size L/2; see 
[22 ] for larger sizes) are shown for conditions I and II at times 
as indicated by arrows. Times tj in the main text appear as 
diamonds on the horizontal axis. All units are arbitrary. 



term becomes effectively positive (and therefore stable) 
at sufficiently large space and time scales [16|, l25| . This 
effect is induced by the nonlinearity and allows for its 
eventual control of the scaling behavior at the stationary 
state, that is in the Kardar-Parisi-Zhang (KPZ) univer- 
sality class [13]. Analogous behavior has been recently 
shown to occur in the 2D isotropic nKS case [175 . Mor- 
phologically, the linear pattern forming instability occur- 
ring at short times is followed by a disordered height 
morphology showing kinetic roughening properties of the 
KPZ class asymptotically. 

In the case of the aKS equation (J2j) with small ax, 
this renormalization can be expected in the x direction, 
resulting in a stabilizing value v^ > 0. In the y direction, 
however, the nonlinearity is so weak that the transition 
to the nonlinear regime is strongly delayed. Therefore, 
at times t c:^ t2 when u^ has already renormalized to a 
positive value, the corresponding Uy coefficient has not 
yet, and remains negative. Then, a new linear instability 
in the y direction causes the formation of a ripple pattern 
that appears rotated by 90° with respect to the early time 
pattern. These ripples then grow exponentially in time 
until the values of the slopes along the y direction become 
so large that the corresponding nonlinearity in Eq. (|2]) 
takes over and the new ripple amplitude is stabilized. 
One can estimate the value of the time ts at which this 
happens to be [23,] 



ts (X a^^ln(a^,/aA) . 



(3) 



This expression arises from the behavior of the surface 
roughness following the linearized equation (J2]) up to 
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FIG. 2: Transition time ts vs (a) a a for a^ — 0.1 and (b) 
au for ax = 10~^. Solid lines represent fits according to Eq. 
([3]) with a prefactor as a fitting parameter. The inset in (b) 
shows the wavelength iy of the rotated pattern as a function 
of azy, as determined from the simulations for ax = 10~^ (blue 
squares; error bars are smaller than the symbol size) and cal- 
culated within a linear approximation (red line) . All units are 
arbitrary. 



times t < ts as W{ts) ^ ex.-p{aiyts/iy)^ the assump- 
tion of a similar scaling for surface height and rough- 
ness, h ^ W^ and a balance between the linear and 
non-linear terms in Eq. (J2j) precisely at time ^3, namely 
ajydyh ^ ax{dyh)^, implying h ~ a^jax. The "linear" 
wavelength of the rotated pattern should moreover be 
given by ^2^ = 27r(2/a^)i/^ 

Quantitative comparisons between the simulated and 
the analytical dependence of ^3 on a\ and a^ are shown 
in Figs. [2fa) and[2jb), respectively. Also, iy is shown 
in the inset of Fig. [2fb) as a function of ay as obtained 
from simulations, and as calculated from the linear ap- 
proximation. In all cases, values from simulations agree 
well with analytical estimates, supporting our interpreta- 
tion of the dynamic morphological transition within the 
framework of parameter renormalization. 

Further progress is possible by studying Eq. ([T]) within 
a one loop Dynamical Renormalization Group (DRG) 
approach. After its application to fluctuating hydro- 
dynamics [26|, this method has recently shown a large 
explanatory power in related contexts, like a multiscale 
descriptions of fluctuating interfaces [23 or morpholog- 
ical instabilities mediated by non-local interactions |28j . 
Following the standard approach [26], we arrive at the 
following KG parameter flow [22], 



dr 
~dl 



v.. 



TV, 



dl 



53 



(4) 



where a coarse-graining of the height and noise fields has 
been performed in a fast mode shell in wave- vector space 
k with k e [A(l — d/). A], where A = 1 is a lattice cut- 
off. Here, r = i^x/^yj 9 = ^l^/i'^'^^y)^ with D being 
the noise variance. Eqs. (|4]) generalize the DRG analysis 
of the anisotropic KPZ equation in [10|. The functions 
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FIG. 3: DRG flow of the couplings r and g. The numerical 
integration of equations Q has been carried out using Ux — 
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of Xy and D as shown in the legends. Upper panel: fixed 
D = 10~^ and changing Ay. Note that g renormalizes even 
for Xy = 0. Markers display values of {r{l),g{l)) separated by 
A/ as shown in the legend, while arrows show the direction of 
the RG flow. Lower panel: for fixed A^ = 0.025, and different 
values of D with A/ = 1/D. All units are arbitrary. 



Tiy. originate from propagator renormalization, whereas 
^ arises in noise variance renormalization [22|. The flow 
(|4|) has been derived within the further assumption that, 
as expected [16, 25], the parameters JCj become enslaved 
[2^ to the slower parameter Uy in Eq. ([T]). In Fig. [3] we 
show results from a numerical integration of @. Starting 
out from conditions for which Ux and Uy are both nega- 
tive, parameter r is seen to cross the zero value meaning 
that, at the corresponding scale /, Vxil) has become pos- 
itive and stabilizing, while Tyy{l) remains negative. After 
further coarse-graining, g decreases to —00, signalling 
renormalization of Uy towards positive values which in 
turn requires crossing Uy = dit the appropriate (large) 
scale. Note that, the larger Xy is, the faster renormaliza- 
tion of g takes place, as denoted by the relative spacing 
among points evaluated at equally spaced positions nAl 
along the corresponding flow trajectories. Once both Uj 
coefficients have renormalized to positive values, by anal- 
ogy with the ID and the isotropic 2D cases for the nKS 
equation one expects the system to enter the anisotropic 
KPZ regime [10]. Within this picture, a stationary state 
is expected at long times, that shows kinetic roughen- 
ing. As in the noisy KS case [15, 16], this state is more 
efficiently reached when the noise variance is larger, as 
seen in the lower panel of Fig. [3l where g renormalizes 
faster for increasing D values. Indeed, the slow growth 
of W{t) at long times for condition II in Fig. [1] signals 
the stabilization of the second ripple structure by non- 
linear effects. Note, this is also the case for condition I 
for which no second ripple structure exists. Thus, bare 
values of the nonlinearities that are comparable to each 



other (condition I) lead to faster renormalization of Ux,y 
to stable positive values and to a rough, disordered sta- 
tionary state, as implied by the upper panel in Fig. [H 

The present dynamic morphological transition induced 
by renormalization effects can be intuitively understood 
by an argument that employs knowledge of a particu- 
lar solution of Eq. ([!]) for the case in which A^y = 0. In 
this extreme limit one can assume [11] that h{x^y^t) = 
H{y, t) does not depend explicitly on x so that d^H = 0, 
even for non-zero u^ and A^^. Substituting this Ansatz 
into Eq. ([T]), we get dtH = i^ydyH — JCydyH + r/, that is 
a linearized ID noisy KS equation for H because Vy < 0. 
Therefore H{y,t) is a ripple structure with crests along 
the X direction. Moreover, Hnear stability analysis leads 
to a dependence of the wavelength of this ripple struc- 
ture precisely as determined for iy in Fig. [2] once rescaled 
coordinates are employed as in Eq. ([2]). Note, how- 
ever, that H{y^ t) is not a solution of Eq. ([1]) for con- 
ditions II. Still, starting out from a flat initial surface, 
the term JVy = Xy{dyh)'^ stays negligible until t '^ t2 
at which the full solution h{x^y,t) « H{y,t). Since the 
latter solution is itself morphologically unstable, large 
values oi dyh ^ dyH build up that make the term My 
no longer negligible. This introduces significant differ- 
ences between h and H after t > t^. A limitation of this 
argument is its neglect of dynamics up to t = ^2 . At ear- 
lier times, a ripple structure exists in the perpendicular 
direction, for which it is ^/-derivatives, rather than x- 
derivatives, that are small (see e.g. morphology Ii, IIi in 
Fig. [T|). For solution H{y^t) to become dynamically rel- 
evant, the initial ripple structure need to be washed out 
by fluctuations and nonlinearity as in condition I, which 
requires the renormalization process discussed above. 

A dynamical transition reminiscent of the one con- 
sidered here has been observed experimentally in high- 
temperature surface nanopatterning by ion-beam sput- 
tering (IBS) of Si(lll) surfaces [3(3, ol]. Specifically, at 
low fluence (equivalent to time for the fixed flux con- 
ditions employed), a ripple pattern with crests perpen- 
dicular to the direction of the incident ion beam formed 
on the surface, with a wavelength i^ — 300 — 500 nm. 
At intermediate fluence, however, a different ripple pat- 
tern rotated by 90° overlayed the initial one, resulting 
in a pattern of dot-like features. At even higher fluence, 
the initial pattern vanished and only the rotated pat- 
tern, with a significantly larger wavelength £y > 500 nm 
remained. The experimentally observed rotation of the 
ripple pattern by 90° does not agree with the predicted 
angle for cancellation modes, which under these exper- 
imental conditions is expected to be 25° [S^. Hence, 



the observed ripple rotation is not related to the appear- 
ance of cancellation modes and must be of a different 
origin. Given the striking similarities to the transition 
studied in our work, a strong nonlinear anisotropy with 
ttA <C 1 can be assumed for the experimental system. 
Transient morphologies of two-dimensional features sim- 
ilar to those observed in [3(/, 31] can also be achieved in 
the simulations of the aKS equation by tuning the a^ and 
ax coefficients in such a way so that the growth of the 
rotated ripples sets in before the initial pattern has fully 
vanished [22|. Although the physical picture leading to 
the aKS equation as a physical model for IBS has been 
recently contested (see a review in [32]), the high tem- 
perature condition employed in these experiments can be 
expected to enhance surface transport. In such a case, an 
aKS-type equation holds with modified coefficients [32[, 
which may moreover account for the lack of cancellation 
modes in the experiments as compared with theoretical 
estimates derived from [8, 30, 31]. 

In summary, we have obtained a dynamical transition 
in the evolution of anisotropic patterns that illustrates 
the rich phenomena that occur when taking into account 
the full anisotropy in this class of non-equilibrium sys- 
tems. The transition merges together two apparently 
opposed phenomena, like the selection of a typical wave- 
length and strong morphological fluctuations leading to 
renormalization and scale invariance. Actually, the dom- 
inance of the latter at intermediate times seem to be a 
requirement for the development of the pattern that later 
emerges. A dynamical role seems to be played also by ap- 
proximate solutions of the equation, such as cancellation 
modes. This suggest the interest of exploring such type of 
solutions in other anisotropic, pattern forming systems. 
In its stabilized form, the isotropic KS equation has been 
shown to provide a generic model for parity-symmetric 
systems featuring a bifurcation with a vanishing wave 
number [33] • Thus, we expect the phenomenology of 
the aKS equation to apply quite generically. In partic- 
ular, the transition that we have discussed may offer an 
explanation for the recently observed ripple rotation in 
high-temperature IBS nanopatterning experiments on Si 
surfaces [30] . Further theoretical and experimental work 
is needed in order to elucidate the degree to which this 
is actually the case, and the appearance of related phe- 
nomena in other pattern forming systems. 
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